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ABSTRACT 

The existence of the Oort Comet Cloud, the Kuiper Belt, and plausible inefficiencies in 
planetary core formation, all suggest that there was once a residual planetesimal disk of mass 
~ 10-100 Mg in the vicinity of the giant planets following their formation. Since removal of this 
disk requires an exchange of orbital energy and angular momentum with the planets, significant 
planetary migration can ensue. The planet migration phenomenon is examined numerically by 
evolving the orbits of the giant planets while they are embedded in a planetesimal disk having 
a mass of Md = 10 to 200 M^. We find that Saturn, Uranus, and Neptune evolve radially 
outwards as they scatter the planetesimals, while Jupiter's orbit shrinks as it ejects mass. 
Higher-mass disks result in more rapid and extensive planet migration. If orbit expansion and 
resonance trapping by Neptune is invoked to explain the eccentricities of Pluto and its cohort 
of Kuiper Belt Objects at Neptune's 3:2 mean-motion resonance, then our simulations suggest 
that a disk mass of order Mo ^ 50 M® is required to expand Neptune's orbit by Aa ^ 7 AU 
in order to pump up Plutino eccentricities to e ~ 0.3. Such planet migration implies that the 
initial Solar System was more compact in the past, with the Jupiter-Neptune separation having 
been smaller by about 30%. 

We discuss the fate of the remnants of the primordial planetesimal disk. We point out that 
most of the planetesimal disk beyond Neptune's 2:1 resonance should reside in nearly circular, 
low-inclination orbits, unless there are (or were) additional unseen, distant perturbers. The 
planetesimal disk is also the source of the Oort Cloud of comets. Using the results of our 
simulations together with a simple treatment of Oort Cloud dynamics, we estimate that ^ 12 
M0 of disk material was initially deposited in the Oort Cloud, of which ^ 4 M© will persist over 
the age of the Solar System. The majority of these comets originated from the Saturn-Neptune 
region of the solar nebula. 

Subject headings: Solar System: formation — Kuiper Belt, Oort cloud 



Introduction 



It is generally accepted that a final distinct stage in the formation of our planetary system consisted 
of the clearing of a residual planetesimal population by the gravitational perturbations of fully-formed 
giant planets. The formation of the Oort Cloud, which is a spherical swarm of comets orbiting the Sun 
at distances of ~ 10^ *° ^ AU, is thought to be a product of this stage of Solar System formation. Mass 
estimates for the Oort Cloud are in the range ^ 10 to lOOAfe ( Weissman 1996 ). Another remnant of 
the planetesimal disk is the recently discovered population of the Kuiper Belt Objects (KBOs) that orbit 
beyond Neptime. Although the mass of the observable portion of the Kuiper Belt is only ~ 0.26 M® 
( Jewitt, Luu. & Trujillo 1998), planetesimal accretion models require an initial mass that is of order ~ 35 
M® for the assembly of Pluto and the R 



100 km-sized KBOs in the 30-50 AU zone (Stern fc Colwell 



1997, Kenyon fc Luu 1998). Noting also that the solid cores of the giant planets are of order ~ 10 Mg 



and that core formation is not likely to have been 100% efficient, it is quite plausible that there was still 
a residual planetesimal disk of mass ^ 10 — 100 M® in the vicinity of the giant planets after they formed. 
The eventual removal of this mass via gravitational scattering by the giant planets could have caused 
significant evolution of the planetary orbits, such that the presently observed orbital configuration of the 
Jovian planets is considerably altered from that which was obtained soon after their formation. 

Not only is an early epoch of planet migration plausible based upon formation considerations, but 
recent advances in our knowledge of the outer Solar System provide new motivation for studying this 
process. The determination of reasonably reliable orbits for several dozen KBOs has revealed that their 
distribution is quite non-uniform: there is a near-complete dearth of KBOs having semimajor axes a 
interior to Neptune's 3:2 mean-motion resonance at 39.4 AU, and there is a prominent concentration of 
objects at the 3:2 resonance with moderately large eccentricities e ^ 0.1 — 0.35 ( Jewitt, Luu, fc Trujillo 
1998| ). Several other KBOs orbit at the 4:3, 5:3, and at the 2:1 resonances, all with moderately large 



eccentricities. There is also a non-resonant KBO population beyond the 3:2 resonance which extends out to 
47 AU; these objects exist in more circular orbits with e ^ 0.1. This peculiar orbital distribution supports 
the hypothesis that Neptune's orbit migrated radially outwards, sweeping the primordial Kuiper Belt with 
that planet's mean-motion resonances, and capturing Pluto as well as a cohort of KBOs at those resonances 
( Malhotra 1993 , Malhotra 1995| ). Both the semimajor axes and the eccentricities of captured bodies would 
have grown concurrently with the planet's orbital expansion. The inclinations of KBOs can also become 
excited when a vertical secular resonance sweeps past ( Malhotra 199q ) . Analysis of the resonance-sweeping 
mechanism shows that if planet migration were responsible for the eccentric and inclined orbits of Pluto 
and the other KBOs at the 3:2 resonance, then Neptune's orbit must have expanded by 5 ^ Aa <, 10 
AU on a timescale of order -r,„ -- lO'^ years or longer ( [Malhotra 1993| , [Malhotra 1998] ). [Malhotra (1997)| 
estimated that the gravitational clearing of a residual planetesimal disk having a mass of about 35 M® 
distributed in the vicinity of Uranus and Neptune would expand Neptune's orbit by Aa '^ 5 AU. The 
planet-migration/resonance-sweeping hypothesis not only accounts for the abundance of eccentric, inclined 
KBOs locked in orbital resonance with Neptune, but it also accounts for the lack of low-eccentricity orbits 
in the 3:2, the depleted region interior to the 3:2 as well as the excited e ~ 0.1 state of the non-resonant 
KBOs. It is thus possible that the characteristics of the orbital migration history of Neptune (and, by 
extension, the other giant planets) is preserved in the details of the KBO orbital distribution. 



Orbital migration was first noted in the giant planet accretion models of Fernandez fc Ip (1984) which 
showed that the orbits of a growing proto- Uranus and Neptune could shift by ~ 5-10 AU while embedded 
in a planetesimal disk having a mass of order ~ 100 M® (Fernandez fc Ip 1984, Fernandez fc Ip 1996). 



Fernandez fc Ip (1984) argue that Uranus and Neptune would preferentially gain angular momentum 



(and thus expand their orbits) by scattering planetesimals and lowering their periheha down to Jupiter 
and Saturn. Jupiter, being an effective planetesimal ejector, would shrink its orbit due to the resulting 
energy and angular momentum losses. However it should be noted that this argument, while plausible, 
is not fully supported by their numerical simulations. Those simulations show net orbital expansion by 
Uranus and Neptune in some instances and orbital decay in others. Individual runs show considerable 
to-and-fro motion by Neptune. Such evolution would only stir up the planetesimals rather than trap them 
at resonances. We note that these simulations were not fully self-consistent since gravitational interactions 
were considered only between planets and planetesimals that were on crossing orbits; long range forces 
between the planets and the planetesimals were neglected, as were the mutual planet-planet perturbations. 
Also, the giant planets' gravitational cross-sections were artificially increased in order to speed up the 
system's evolution. This approximate treatment of planet-planetesimal dynamics was necessitated partly 
due to the limited computer resources available at the time. A more realistic treatment of disk clearing and 
planet-migration is clearly warranted. 

The present work revisits the planet-migration scenario using more sophisticated numerical techniques 
to evolve a system of giant planets that is embedded in a disk composed of many low-mass particles. In 
Section 0, we describe the model and test results in detail. Primary results on the giant planets' orbital 
evolution with disks having different initial masses are given in Section pi Since the planetesimal disk is 
the source of Oort Cloud, we discuss the implications of our results for the formation of the Oort Cloud 
in Section ^. Particles that manage to avoid ejection will reside in what are sometimes called the classical 
and the scattered disks, which are characterized in Section |[ The role of unmodeled effects such as disk 
self-gravity and the role of spiral density waves are discussed in Section |6|, and conclusions are given in 
Section 0. 



2. Numerical Method 

Our model consists of the Sun, the four giant planets Jupiter, Saturn, Uranus, Neptune, and a 
population of numerous (N) low-mass particles initially distributed in a disk. Ideally, one would like to 
have N ~ 10"'^° to follow the evolution of this system with a completely self-consistent calculation of the 
mutual gravitational perturbations and also include the external forcing due to the galactic tide. However 
such a calculation exceeds current computational resources. Instead, as we describe below, we use a fast 
orbit integrator together with several simplifications to obtain an approximation to the ideal within a 
reasonable amount of computing time. 

In our simulations, all bodies are treated as point masses, thus the possibility of collisions or accretion 
is neglected. This approximation is justified during the late stages of giant planet formation when the 
frequency of planet-particle collisions is smaller than the scattering frequency by a factor of the planet's 
physical radius/Hill radius squared, which ranges from 10~^ for Jupiter to 10^* for Neptune. 

The model includes the mutual gravitational forces exerted between the planets as well as the forces 
between the planets and the low-mass particles, but the forces amongst the particles themselves are 
neglected. This approximation is employed to reduce the computational expense, and it also eliminates the 
unphysical self-stirring that would otherwise occur in the model disks considered below that have masses of 
Md = 10 to 200 Mg distributed among just N — 1000 particles. Had the particle-particle interactions been 
included, the relatively small number of N massive disk particles would quickly stir themselves up so much 
that the resulting system would hardly resemble a planet-forming particle disk. Although this difficulty is 



avoided by neglecting the particle-particle forces, this approximation also precludes the possibility that the 
planets might generate spiral density waves in the disk (see Section y) . 

Dynamical models have shown that particles scattered into wide orbits with semimajor axes a ^ 3000 



AU become decoupled from the planets due to the galactic tide (Duncan, Quinn, & Tremaine 1987; see also 
Section ||) . These distant bodies are deposited in the Oort Cloud, which is effected here by simply removing 
particles from the system when a > 3000 AU. 

A second-order mixed variable symplectic (MVS) mapping is used to rapidly advance the heliocentric 
positions and velocities of the planets and low-mass particles as they interact in the Sun's gravitational 



field. Our implementation follows the algorithm of Wisdom & Holman (1991) with the improvements of 




Saha fc Tremaine (1992) . A fixed step size of At = 0.4 years is used, which is sufficiently short to resolve 
the orbits of particles that evolve down to perihelia of q ~ 3.5 AU as well as most encounters between 
particles and planets. For the latter, the dynamical timescale of a typical scattering event isn 

(1) 

where Mp is the scattering planet's mass and q is the particle's closest approach distance. Evaluating T* 
with q set to each planet's Hill sphere radius shows that T* ~ 4.1 years at Jupiter, 10 years at Saturn, 29 
years at Uranus, and 57 years at Neptune. Thus the integration timestep of At = 0.4 years is at worst 1/10 
of the dynamical time for orbits that graze Jupiter's Hill sphere. However an MVS algorithm that employs a 
fixed step size will fail to resolve very close planet-particle encounters that penetrate well within a planet's 
Hill sphere, and will also fail to correctly evolve a very eccentric orbit when the perihelion passage is not 
well time-sampled ( Ranch fc Holman 1999| ). These difficulties are mitigated by the using the following 
approximat ions . 

(i) During a close encounter between a particle and a planet, two-body trajectories are adopted for 
their relative motion whenever a particle passes sufficiently close to a planet. Experimentation shows that a 
"sufficiently close" encounter is one that changes the particle's fractional distance from the planet by more 
than 50%, or varies its planet-centered angular coordinate by more than 90°, during the time step At. For 
particles that are initially in low-eccentricity heliocentric orbits, the two-body approximation is triggered 
only when a particle approaches within about 10% of the planet's Hill sphere. For more distant encounters, 
particle trajectories are in fact evolved with greater accuracy using the standard MVS mapping. 

(ii.) When a particle evolves into an orbit of sufficiently high eccentricity, its motion during perihelion 
passage may be poorly time-sampled and the MVS mapping can produce unphysical evolution. For the 
step size employed here, significant errors accrue in orbits having perihelia g ^ 3.5 AUFl. Proper treatment 
of such orbits is not of minor import since approximately 30% of the disk particles cycle through the inner 
Solar System with heliocentric distances of r < 3.5 AU. A simple, but impractical solution, is to use a step 
size small enough to resolve all perihelion passages. However a more practical approach consists of turning 
off all planetary perturbations while a particle travels interior to a heliocentric distance of r < 3.5 AU; 



^This encounter time is defined as the time to change the particle's true anomaly from — 7r/2 to 7r/2 for a planet-centered 
parabolic orbit having a pericenter distance q from the planet. 



It should be noted that the symplectic integrators of Levison & Duncan (1994) and Duncan, Levison, & Lee (1998) 



do not 

suffer this instability when massless test particles achieve low perihelia orbits (H. Levison, personal communication). However 
the intervention described above is still required since there is no fully symplectic algorithm that can compute low-perihelia 
trajectories when the particles carry mass. 



this results in a piecewise Keplerian trajectory that is approximately correct. With this approximation, 
the particle's perturbations upon the planets are still not fully time-sampled, which causes slow drifts in 
the system's total energy and angular momentum. Numerical experimentation has shown that the most 
economical and robust solution to this problem requires turning off the particle's perturbations upon the 
planets while its perihelion distance is q < 3.5 AU. 

The remainder of this section describes a few critical tests of the algorithm in order to demonstrate 
that the approximations made here do not introduce any artificial orbital evolution of the planets and the 
particles that is of any significance. 

The quality of the orbit integrations may be illustrated with the restricted three-body problem that 
consists of a massless particle perturbed by a planet on a circular orbit about the Sun. For this system the 
particle's Jacobi integral J = E — L^^p is conserved, where E is the particle's energy, L^ is its component 
of angular momentum perpendicular to the planet's orbit, and fip is the planet's mean motion. Figure n^ 
shows the mean fractional variations AJ/J that result when integrating 100 particles for 10* years with 
a Jupiter-mass planet at a^ = 5 AU. Squares are shown for those particles that start in rather circular 
orbits near and beyond the planet; these experience a AJ/J < 10~^ that decreases with distance provided 
they have not already made a close approach to Jupiter. Filled squares indicate particles that approach 
closer than a Hill distance to Jupiter; these experience larger variations in J due to the close-encounter 
approximation used here. The remaining particles (crosses and circles) are all on Jupiter-crossing orbits 
with their perihelia ranging over 0.05 < q < 6 AU. Crosses are shown for particles that do not encounter 
Jupiter, whereas filled circles represent particles that do encounter the planet. The crosses show that 
AJ/J grows inwards of Jupiter's orbit due to the poor time-sampHng of the particles' perihelion passages. 
However the error AJ/J levels off where the integrator turns off the planet's perturbation and the particle 
motions are temporarily Keplerian. Failure to implement this procedure would have instead yielded 
disastrous results inwards of about g ~ 1 AU. Shrinking the step size At would of course reduce the error 
in J but at the expense of precious computer cycles. 

An alternate scheme that avoids the growth of numerical errors among low-g particles is to simply 
remove them when they drop down to a heliocentric distance of r < 3.5 AU. Although this procedure 
inhibits error growth, it is less than desirable since it introduces an unphysical (and very hungry) mass sink 
to the system. When employing this procedure to the simulations of planet-migration described in Section 
S (which employ the Keplerian approximation whenever r < 3.5 AU), we find that the orbital migration of 
Saturn, Uranus, and Neptune proceeds radially outwards at similar rates. However the removal of r < 3.5 
AU particles results in an artificial sunward mass-flux that causes Jupiter's orbit to expand rather than 
shrink. Nonetheless, the fact that the orbital evolution of the outer giant-planets' is rather insensitive to 
the treatment of low-q particles adds confidence that the Keplerian approximation employed here is not 
driving the planet-migration reported in Section y. However enforcing the particles' Keplerian motion 
inside of r < 3.5 AU is still needed in order to accurately model Jupiter's orbital history. 

In the absence of any close encounters between particles and planets, the MVS integrator preserves the 
system's integrals in the expected manner. However Fig. [l| shows that the close encounter algorithm used 
here will not preserve a particle's J to better than 1 part in 10"^ should it pass nearer than a Hill distance of 
a planet. This numerical error results in an unphysical diffusion of particle trajectories. This is a concern 
since the planet migration phenomenon is sensitive to the flux of particles that are scattered and exchanged 
amongst the planets (Fernandez & Ip 1984, Malhotra 1993| ). 



To judge the effects of numerical diffusion, the restricted three-body problem is again integrated for a 



system of 100 massless particles that orbit in the vicinity of a Neptune-mass planet in a circular orbit at 
Up = 30 AU. The particles' initial semi-major axes range over 28.7 < a < 31.5 AU and have small initial 
eccentricities e and inclinations i. For this system, each particle's dimensionless Jacobi integral, 



— {l-e^)cos^ + Oi^l), (2) 

Up 

is a conserved quantity. Here, a particle's semimajor axis a is in units of Up and /i — 5.15 x 10^"' is Neptune's 
mass in solar units. In this test run, the particles initially have a mean Jacobi integral (J) ~ 2.998 with 
a dispersion aj — 5 x 10^^. To obtain a physically meaningful measure of the numerical diffusion that 
occurs in the model, define a threshold value Ju — 2.989 which corresponds to an orbit having a perihelion 
at 20 AU (near Uranus' orbit) and an aphelion at 30 AU (near Neptune's orbit) so that a/up = 5/6 and 
e — 1/5. Although this orbit is dynamically forbidden to the test particles, numerical diffusion might 
allow particles to cross Uranus' orbit which would result in an unphysical exchange of particles between 
Neptune and Uranus. The system is integrated for 5 x 10^ years, and at the end of the run 11 particles 
have diffused into forbidden Uranus-crossing orbits that have J < Jjj- The timescale to diffuse into crossing 
orbits, td, is obtained from Fig. which shows the particles' average J and their dispersion <jj versus time. 
The characteristic diffusion timescale for particles to diffuse across the Jacobi 'gap', i.e., the time when 
I J — <Tj\ < Ju, is td ~ 2.5 X 10^ years. 

However it should be noted that a particle's J is not conserved in a multi-planet system. For example, 
adding a Jupiter at 5 AU to the above simulation will drive these particles into Uranus-crossing orbits 
having J < J[/ at a rate that is about 10 times faster. Since the numerical diffusion rate is considerably 
slower than the dynamical diffusion that occurs in a multi-planet system, we conclude that the transport 
of particles between the planets due to numerical diffusion is not significant in these simulations. 

As already noted, it is the close encounters between the planets and the planetesimals, and their 
concomitant exchange of angular momentum, which drive the planet migration process. In our simulations, 
the planet-particle relative motion during very close encounters is not computed exactly, so any small 
error in the angular momentum exchange is also reflected in the recoiling planet's orbit. The final test 
discussed below verifies that the errors in this angular momentum exchange are in fact too small to drive 
the planet-migration described in Section 0. 

The motion of Neptune as well as several hundred particles of infinitesimal mass m are integrated at 
the usual step size At — 0.4 years for 200 years. The same experiment is then repeated using a step size 50 
times smaller. Since integration errors decrease with step size, this latter run may be regarded as a much 
more exact representation of the particle trajectories. Upon differencing the two runs, the error in each 
particle's velocity relative to the planet, (5V, is calculated as the particle exits the planet's Hill sphere (see 
Fig. pK). The Figure shows that most encounters occur at the periphery of the planet's Hill sphere which 
result in rather small relative velocity errors. However the very close encounters having periapse q ^ O.IRh, 
which account for 6% of these encounters, result in sizable velocity errors that are of order 6V/V ~ 0.1. 
The error in the specific angular momentum £p exchanged between the planet and a scattered particle is 
Sip — —mVp X (5V/(Afp + m), where r,, is the planet's heliocentric coordinate and Mp is its mass. Figure ^b 
shows a histogram of the z-component of the errors Sip which are distributed about zero. The net specific 
angular momentum exchanged between the planet and the particle swarm is the sum L = | ^€p|, and its 

root-mean-square error is 6L — ^/X^'^^p- For the encounters shown in Fig. y, the fractional error in the 
angular momentum exchanged is only 5L/ L — 0.016. The above procedure is also repeated for encounters 
at Jupiter, which are less well time-resolved, and yields a larger fractional error 5L/L = 0.12. It should 



be noted that these fractional errors decrease for higher relative velocities, which reflects the fact that the 
two-body approximation, which neglects the Sun's gravity, becomes more accurate at faster encounters. 
Thus the error in the total angular momentum exchanged between the planet and neighboring particles 
will steadily decrease as the planet heats up the particle disk. Since these fractional errors are small, we 
conclude that the close encounter approximation used here does not drive the planet migration described 
below. 



3. Simulations of Planet Migration 

3.1. Initial conditions 

In the simulations reported here, we adopted initial planet orbits similar to those used in previous 
investigations of planet migration, in which the semimajor axes of Jupiter, Saturn, Uranus, and Neptune 



are displaced from their present orbits by respective amounts Aa = +0.2, -0.8, -3.0, and -7.0 AU (Malhotra 



1995). The planets are assumed to have their present masses. The initial planetesimal disk is composed of 
1000 equal-mass particles distributed in orbits of 10 < a < 50 AU such that the disk's inner edge lies just 
exterior to Saturn and the outer edge lies just beyond the present location of Neptune's 2:1 resonance. The 
disk surface density a varies as a~^ , a{a) = 4.0 x 10~^{Md/A[J )(1 AU/a), where Md is the total disk 
mass. Thus, approximately two-thirds of the disk starts exterior to Neptune's initial orbit. Four separate 
simulations are presented below in which the disk's initial mass is Mjj = 10, 50, 100, and 200 M^, and the 
individual disk particles have masses m — 0.01, 0.05, 0.1 and 0.2 M®. An additional 50 massless particles 
are also distributed between 50 < a < 100 AU in order to assess the degree of perturbation of a hypothetical 
part of the Kuiper Belt extending well past that which is presently observable. 

It should be noted that accretion models advocate an initial disk containing several tens of Earth-masses 



in the 30-50 AU zone in order to form Pluto and QBi-type Kuiper Belt Objects (Stern & Colwell 1997, 



Kenyon & Luu 1998), a scenario that is bracketed by the models explored here. 



For those particles that initially lie far from any planets, initial eccentricities of e^ = 0.01 and 
inclinations i^ = 6^/2 are adopted. However particles having a semimajor axis near a planet will already 
have experienced a history of stirring that results in particle dispersion velocities Vd of the order of the 



planet's escape velocity at its Hill sphere radius Rh ( Ida fc Makino 199S ), where 

Vd -- y/2GM/RH and Rh = (M/3M0)^/^a, (3) 

G is the gravitation constant, M the planet's mass, and Mq the solar mass. Assuming 
the particles' inclinations are half their eccentricities, the particles' dispersion velocity is 



Vd — y/i^d + i'd)GMQ/a ^ Cdy/bGAlQ/Aa which corresponds to an eccentricity of ed ^ y/24j5{RH / a) . 
With Rh /a ~ 0.025 for both Uranus and Neptune, the particles starting in each planet's heated zone have 
initial e^ — 0.06 and id = 1.7°. The adopted half-width of each planet's heated zone is simply each planet's 



feeding zone' of Aa = 2^Rh (e.g. |[da fc Makino 1993| ) such that Aa = 1.4 AU for Uranus and 2.0 AU for 



Neptune when evaluated at their initial heliocentric distances. 

For a planet that is embedded in a swarm of identical particles of mass m, dynamical friction will tend 
to seek an equipartition of energies in the system's epicyclic motions. Again assuming a planet's inclination 
obeys ip = ep/2, its initial eccentricity is then Cp = edym/M, where e^ are the eccentricities of particles in 
the heated zone. Since Jupiter and Saturn start interior to the particle disk, initially circular and coplanar 
orbits are adopted for these two planets. 



3.2. Results 

Figure shows the orbital histories of the four giant planets as they scatter the surrounding disk 
particles in each of the four simulations. Wc find that the lowest-mass disk, M]j = IOMq, yields little 
evolution in the planets' orbits, but the higher-mass disks result in significant radial displacements of 
the planets during the 30 Myr runs. These simulations confirm the expectation that Saturn, Uranus and 
Neptune migrate radially outwards while Jupiter migrates slightly inwards. For a given disk mass, the 
magnitude of radial migration is largest for Neptune and successively less so for the interior planets. As 
might be expected, planetesimal disks of greater mass result in planet migration that is larger in magnitude 
and more rapid, and also more stochastic owing to the individual particles' greater mass. Although the 
planets' orbital eccentricities and inclinations remain small, there is a clear trend towards larger e's and i's 
at larger disk masses. Since their final e's and z's are largely determined by the numbers and masses of disk 
particles used in the simulations, their final state have little relation to the giant planets current e and i 
configurations. It is worth noting that, in the high-mass disk simulations, the planets pass through a few 
mutual low-order mean-motion resonances, but the planets' orbits do not persist in any resonance-locked 
configurations. This is not entirely surprising as the general trend in the orbital evolution is such that 
the planets are driven towards greater mutual orbital separation, which is not conducive to maintaining a 
resonance libration ( Dcrmott, Malhotra fc Murray 198q ). 



Figure H shows the state of the AIjj — 50 Mq system at logarithmic time intervals. Large dots indicate 
the planets' eccentricities and semimajor axes while small dots/crosses denote particles that have/have not 
had a close approach to a planet. Those particles scattered by the planets yet still bound to the Sun tend 
to have perihelia between the orbits of Saturn and Neptune, as indicated by the two curves. It should be 
noted that in all of the simulations reported here, the orbital migration of the planets has not ceased by 
the end of the runs, and that further planet-migration will continue (albeit more slowly) past 30 Myr. 
The simulation of the Md = 50 M^ disk is extended to 50 Myr in Fig. H, which shows Neptune slowly 
expanding its orbit out to 30 AU. These results suggest that in order to actually 'park' Neptune at a = 30 
AU requires an adjustment of parameters towards a slightly lower disk mass, and/or a steeper gradient in 
the disk's surface density profile, and/or an outer disk edge that lies closer than 50 AU. 



3.3. Discussion 

The disk particles used in all of these simulations are sufficiently massive that their perturbations 
upon Neptune result in non-adiabatic expansion of that planet's orbit. However it is evident in Fig. ^ that 
planet migration is smoother when the disk particles have lower mass. Thus a more realistic simulation 
employing larger numbers of particles will better resolve the disk-planet perturbations. We expect that such 
simulations will exhibit orbit migration which proceeds more smoothly. Nearly adiabatic orbital migration 
is in fact required if Neptune is to efficiently capture particles at its exterior mean-motion resonances. Since 
smooth outward expansion was not realized in these simulations, Neptune did not capture any particles 
at its mean-motion resonances. Thus the simulated disk's end state cannot be directly compared to the 
delicate resonant structure that is observed in the Kuiper Belt. Nevertheless, these simulations do provide 
useful constraints on the likely mass of the initial debris disk that may have been present during an early 
epoch of planet migration. 

Figure reports the radial displacement Aa experienced by each planet after i = 3 x 10^ years as a 
function of the initial disk mass Mo- As noted above, Jupiter migrates sunwards while the other planets 



10 



migrate outwards in each simulation. The time-averaged torque Tq — AL/t that drives Neptune by 
increasing its angular momentum by AL during the run-time t is shown in Figure ^. This torque, as well as 
the displacements Aa of Fig. ^ should be regarded as upper limits since resonance trapping did not occur 
in these simulations. Had resonance capture been realized here, an opposing torque on the planet would 
have developed since the planet must transfer angular momentum to expand the orbits of particles trapped 
at resonances. 

The back-torque due to a ring of mass m that is captured at an exterior j + 1 '■ j mean-motion 



resonance can be calculated from the rate of change of its angular momentum, L,„ = my^GMQa{l — e^). 
Since the resonant torque from the planet expands the ring's semimajor axis at the rate a/a = dp/ap and 



also pumps up eccentricities at the rate de^ /dt — (dp/ap)/ {j + 1) (Malhotra 1993), the rate of change of Lm 

is 

dim m 

—— ~ a T (4) 

at rup 

where a = ap/a = {j /{j + 1))^^'^ is the ratio of the semimajor axes of the planet and the resonance site, 
rUp is the planet's mass, and T — rup^j GM^apdp/^ap is the net torque on the planet, which simplifies 
to T == To ^ dLra/dt ~ ro/(l + am/mp). Thus the back-torque from the resonance captured particles 
will significantly slow the planet's orbital expansion if the mass trapped at resonance is in excess of the 
planet's mass: m ^ rrip/a. In a generic system, the strongest exterior mean-motion resonance is the 2:1 
(j = 1, a = 0.63), which can be expected to capture the most planetesimal mass. If the capture efficiency is 
£, then we can estimate that planet migration will slow if 

£MD,ext ^ l-Gnip (5) 

where M^^ext is the mass of the disk exterior to the planet's initial orbit. 

If resonance capture is rather inefficient, then a high-mass disk is simply fuel for planet-migration. 
But if resonance trapping is effective, say, e ~ 50%, then the planet's migration will slow after the planet 
has captured a mass m comparable to its own at its exterior resonances. In the disk models considered here 
(with Mrt.ext ~ '^/sMjj), this would occur only in a high-mass disk with Af^i > 80 M0. 

If migration by Neptune across a distance Aa ~ 5-10 AU is to explain the sculpted appearance of 
the Kuiper Belt, then Figs. ^ and suggest that the initial planetesimal disk must have had a mass Ado 
in excess of 10 M^ (since a lower-mass disk results in insufficient planet migration that proceeds too 
slowly) but likely less than ^ 100 M^ (since such a high-mass disk would likely produce additional giant 
planets). A more precise disk-mass estimate requires detailed knowledge of the disk's radial extent and the 
particulars of the disk surface density variations a (a). 



4. The Oort Cloud Mass 

The observed flux of long-period comets provides a constraint on the present mass of the Oort Cloud. 
In this section, we use the results of our numerical simulations to estimate the mass of its progenitor, the 
primordial planetesimal disk in the outer Solar System. 

After adjustment for the efficiency of comet detections it is estimated that ^ 63 long-period 



comets/year pass through the inner Solar System with perihelia g < 4 AU (Everhart 1967). This flux stems 



entirely from the Oort Cloud; objects originating in the Kuiper Belt evolve instead into the short-period 



Jupiter-family comets (Levison 1996). Dynamical models of the Oort Cloud, when adjusted to match the 
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flux of long-period comets, require a reservoir of Nqc ^ 1-6 x 10^^ bodie^ ( [Hcislcr 1990 ). Multiplying 



by the typical comet mass yields the total mass of the Oort Cloud. Based upon an admittedly uncertain 



relationship between cometary brightnesses and their size, Weissman (1996)| concludes that the mean comet 
mass is of order ~ 10^® gm, indicating an Oort Cloud mass of order Mqc ~ 27 M^. However this mass 
estimate should be regarded as uncertain by at least an order of magnitude since it relies upon a host of 
uncertain quantities such as the efficiency of comet detections, the comet brightness-mass relationship, 
as well as uncertainties in the Oort Cloud perturbations (e.g., the strength of the galactic tidal field, the 
frequency of stellar encounters, etc.). 



The dynamics of the Oort Cloud is succinctly summarized by Duncan, Quinn, & Tremaine (1987). 
Unless a particle is otherwise ejected from the system, planetary perturbations cause its semimajor axis to 
diffuse both inwards and outwards while keeping its perihelion locked in the giant planet region of the Solar 
System. However a more distant particle is susceptible to perturbations by the galactic tide and passing 
stars which cause its perihelion to diffuse on a timescale that varies as tq oc a^^ . Those particles reaching 
3 X 10'^ ^ a <; 2 X 10* AU which have had their perihelia raised well beyond the orbit of Neptune are thus 
decoupled from the planets and are usually identified as inner Oort Cloud comets. However the perihelia 
of more distant bodies diffuse more rapidly, and those with 2 x lO'' <^ a < 10^ AU that reside in the outer 
Oort Cloud are in fact more likely to diffuse back into the inner Solar System and become observable as 
new long-period comets. It is this flux of new comets that provides an important constraint on the amount 
of mass driven from the initial disk that has managed to avoid ejection during the last 4.5 Gyr. 

Estimates of the Oort Cloud mass may be inferred from the data given in Table n^. The quantity f^k is 
the fraction of the dynamically active disk that diffuses into the Oort Cloud at a > 3 x lO'^ AU during each 
simulation. The dynamically active disk is that part of the disk where particles are likely to be perturbed 
into Neptune-crossing orbits over the age of the Solar System. Long-term integrations of Kuiper Belt orbits 
show that particles in the active disk have semimajor axes a <J l.AaN, where qn is Neptune's semimajor 



axis ( Duncan, Levison, fc Budd 1995 ). The mass of the dynamically active disk. Mad (Table |i|), is defined 



here as the total mass of ejected particles plus all survivors having a <J 1.4aAr that are presumably in 
unstable orbits. Also given is fh, which is the fraction of the active disk that has been ejected from the 
system. Fig. O displays fsk and fh versus time for the Mo = 100 M^ run. What is most striking is that 
at the end of all four runs both f^k and fh vary little among the different simulations (Table ||). This 
indicates that the total mass deposited in the Cloud depends only on the mass of the disk that lies within 
the planets' gravitational reach, and is not very sensitive to the orbital histories of the migrating planets. 

Ultimately, all particles starting in the dynamically active disk are either ejected, deposited in the 
Oort Cloud, or in some instances accreted by the planets. Although the latter outcome is not modeled 



here, impacts may be assessed ex post facto using the collision probabilities of Opik (1951) . After summing 
the probability of each particle striking each planet, we find that ~ 15 of the 1000 particles would have 
struck the giant planets during each simulation, with roughly half of these impactors striking Jupiter. Such 
impacts would have contributed no more than ^ 2% to any planet's mass, so the neglect of particle-planet 
collisions is justified. It is also worth noting that a few percent of the disk passes through the terrestrial 
zone, as is indicated by the /go.5 curve of Fig. S which shows the instantaneous disk fraction having 
perihelia inside of 3.5 AU. In these simulations typically ~ 30% of all disk particles have brief episodes with 
q < 3.5 AU. These findings are consistent with earlier studies showing that planetesimals scattered during 



3n 



The Heisler 199C Oort Cloud model shows that a reservoir of 10^^ comets will produce a flux of 0.2 comets/year having 



q < 2 AU. To match the observed flux of comets with g < 4 AU, multiply this population by 63/(2 X 0.2). 
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the epoch of disk-clearing and planet-migration may have contributed significant numbers of impactors 



during the late heavy bombardment of the terrestrial planets ( Wetherill 1975 , Shoemaker fc Wolfe 1984 ). 

Table n^ shows that in the four simulations the active disks are depleted by a factor / = f^k + fh — 50% 
after i = 3 x 10^ years. It should be noted that the family of f{t) curves is described well by a power-law 
/ ex i°-^^. This indicates that a fully evolved system having / ^ 1 requires an integration lasting the 
duration of the disk's dynamical lifetime t^ ~ 1.5 x 10^ years, which is beyond our computational means. 
However it is straightforward to extrapolate the formation of the Oort Cloud from the simulations at hand. 
Fig. y shows that the ratio fak/fh ~ 0.4 remains relatively constant during the bulk of the run. This 
relation permits an extrapolation to a fully evolved state of the system having /g^ + //^ -^ 1, where the 
primes denote final extrapolated values. Assuming f^i^/fl^ = fsk/fh over the age of the solar system, the 
extrapolated fractions become /^ = (1 + fsk/ fh)~^ — 0.73 and /g^, = 1 — /^ ~ 0.27; these values are also 
given in Table H. Since these disk fractions are both number as well as mass fractions, the extrapolated 
planetesimal mass that is initially deposited in the inner Oort Cloud is f^i^Mad- However the galactic 
tide will subsequently strip away comets that diffuse past a ^ 10^ AU, and passing stars will eject others. 
Numerical studies show that these external perturbations acting over the age of the Solar System will 



reduce the Oort Cloud mass to about a third of its initial value (puncan, Quinn, fc Trcmaine 1987). Thus 



the final Oort Cloud mass reported in Table || is Mqc — fsk-^da/S, and is also displayed as the solid curve 
in Fig. |l^ as a function of the initial disk mass Md- 

The nebula origin of the Oort Cloud is given by the solid curve in Fig. pi which shows a histogram 
of the Oort Cloud particles' initial semimajor axes for the AId = 50 M^ disk; after t = 3 x 10^ years, 
the planetary configuration in this system most resembles the Solar System. Since the inner edge of our 
model disks is truncated at 10 AU, any Oort Cloud mass originating in interior orbits is still unaccounted 
for. We estimate this contribution after the fact by evolving a system of four giant planets in their present 
configuration plus an annulus of 153 massless test particles having an initial a ex. a~^ surface number 
density between 4 < a < 10 AU. Their Oort Cloud contribution is given by the dashed curves in Figs. |lO| 



and 11, Evidently, all parts of the giant planet domain contribute mass to the Oort Cloud. These findings 



are in agreement with Weissman fc Levison (1997) who reported that a small fraction of Oort Cloud bodies 
can originate from orbits interior to Jupiter and thus have asteroidal rather than cometary compositions. 

If the initial disk had a mass 10 <; Md <; 100 M^ (Section y|), then Fig. 10 indicates that the resulting 



Oort Cloud mass is 0.5 <^ Mqc ^ H M®. This mass estimate must be qualified for two reasons. For the 
lower-mass disks, Mo < 50 M^, the mass estimate Mqc as given in Fig. |l^ is likely an underestimate, 
because in deriving it we have not accounted for the possibility that Neptune can migrate deeper into the 
disk as the system evolves further on timescales longer than ^ 50 Myr, thus allowing additional material 
to be injected into the Oort Cloud. (This does not affect the higher disk-mass simulations since Neptune's 
gravitational reach has already swept across the entire disk in the higher-mass runs.) For the higher-mass 



disks, Mqc is likely overestimated in Fig. 10 because resonance trapping tends to slow planet-migration 



and may reduce the mass encountering the planets and thus the mass deposited in the Oort Cloud. We note 



that an Oort Cloud having a total mass Mqc ^ 5 Mg and a population Nqc ^ 1-6 x lO^'^ comets (Heisler 



1990| ) suggests that a characteristic Oort Cloud comet has a radius ^ 1 km for a density ^ 0.5 gm/cm' 



3 
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5. The Resonant, Stirred, and Scattered Kuiper Belt Components 

Although the giant planets scatter planetesinials throughout the entire Solar System, a large gap in 
orbital phase space will develop, as illustrated by Fig. |l^. This gap is easily explained via the restricted 
three-body problem, which shows that a planet on a circular orbit will scatter a massless body in a 
manner that preserves the particle's Jacobi integral J (Eq. |^). Particles originating in a cold disk in the 
vicinity of the planet have semimajor axes a "^ ap and J ~ 3. Once scattered, these particles will have 
eccentricities defined by the curve ej=3 (Fig. O), assuming i = 0. It has also been shown that a planet 



of low eccentricity will repeatedly scatter particles along a curve that approximately preserves J ~ 3 (Ida 



fe Makino 1993). Thus when several planets are present, those interior to Neptune can loft particles into 



high-eccentricity Neptune-crossing orbits and ultimately fill the e > ej=3 region of phase space. Those 



particles in the high-eccentricity e ^ ej^a orbits of Fig. 12 are referred to as the scattered disk (e.g.. 



Duncan fc Lcvison 19971) . The region occupied by bodies of lower eccentricity beyond ^ 35 AU is sometimes 
referred to in recent literature as the 'classical' Kuiper Belt; the eccentric, resonant objects (including the 
so-called Plutinos at Neptune's 3:2 resonance) may be considered a distinct dynamical sub-class of the 
latter population. 



Although Fig. ha shows the scattered and classical disks after only 5 x lO"^ years, Duncan fc LevisoE 



1997)1 have integrated test particle orbits for the age of the Solar System. They find that while most 
particles in the scattered disk are removed in less than the age of the Solar System, about 1% of the 
scattered particles survive for longer times (perhaps by acquiring protection via the Kozai mechanism 



(Kozai 1962) or by sticking near mean-motion resonances). In the absence of disk-stirring by any other 
large distant perturbers, we note that the phase-space gap between the classical and scattered disk will 
persist over the age of the Solar System. This gap should become evident as deeper observations begin to 
peer beyond Neptune's 2:1 resonance. 

Until very recently, no KBOs were known to orbit at Neptune's 2:1 resonance, in apparent conflict with 
the prediction of the planet-migration/resonance-sweeping theory that KBO populations at the 2:1 and 3:2 
should be similarly abundant and have comparable eccentricities ( Malhotra 1995| ). As this paper was being 



written, we learned that two KBOs, 1997 SZig and 1996 TRge, have been identified as librating in the 2:1 



Neptune resonance (Marsden 1998|). We note that these two KBOs were formerly identified in Neptune's 5:3 



and 3:2 resonances, and that their orbits were revised to the 2:1 resonance only after observations spanning 
two and three oppositions, respectively. Clearly, orbit-fitting biases and observational incompleteness 
remain in the current census of the Kuiper Belt. A robust test of the planet-migration/resonance-sweeping 
theory requires a larger observational sample of reliable KBO orbits. 

I I 

Morbidelli & Valsecchi (1997) offer an alternative explanation for the structure in the Kuiper Belt. 

They suggest that if a stationary Neptune had scattered a couple of Earth-mass planetesimals (e.g., LNSPs 
= large Neptune-scattered planetesimals) outwards, these massive bodies could have stirred up KBO 
eccentricities to e '--^ 0.2, similar to those observed for KBOs at Neptune's 3:2 resonance. Gravitational 
scattering can indeed insert KBOs into, as well as remove objects from, mean-motion resonances, but 
scattered objects tend to librate at resonance with such large amplitudes that close encounters with 



Neptune become possible and long-term orbital stability is precluded (Levison & Stern 1995). Therefore 
an additional sequence of collisions and/or scattering events is required in order for particles to diffuse to 
stable, low-amplitude librating orbits. Only the fortunate few would survive this process, so the yield of 
KBOs scattered into stable e ~ 0.2 orbits at the 3:2 resonance would be much smaller than that which 
might otherwise be acquired by means of adiabatic orbit expansion and resonance-sweeping by Neptune. 
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Morbidelli fc Valsccchi (1997)| also point to the e ^ 0.1 KBOs that reside between the 3:2 and 2:1 resonances 



as additional evidence for LNSPs. But resonance capture is not entirely efficient, and similar eccentricities 
are also achieved as Neptune's orbit expands and its 2:1 resonance sweeps across the disk and stirs up 
eccentricities (Fig. [l^ ). 

A census of KBOs beyond Neptune's 2:1 resonance (at 48 AU) would permit an evaluation of possible 
stirring by hypothetical LNSPs. If the natal planetesimal disk does extend past 48 AU, then Fig. |2| shows 
that planet migration will produce a stirred zone interior to the 2:1 yet leave the disk exterior to the 2:1 
relatively undisturbed. However, an abundant population of eccentric KBOs beyond the 2:1 would suggest 
a history of additional stirring by other unseen perturbers (though this would not preclude an episode of 
planet-migration) . 



6. Disk Self— Gravity and the Role of Spiral Density Waves 

As noted earlier, our simulations neglect the planetesimal disk's self-gravity in order to inhibit an 
unphysical degree of self-stirring. The consequences of disk self-gravity could be better studied only by 
simulating disks composed of many more lower-mass particles. In this section, we discuss the possible 
consequences of the disk's self-gravity. 

If the local disk mass exceeds the mass of a nearby planet (as is the case for Neptune in most of the 
simulations considered here), it is the disk's gravity that can control the rates at which the perihelia and 



nodes of both the planet and the disk particles precess and/or regress ( |Ward 1981| ) . Secular resonances are 
sites in the disk where a planet's perihelion/node varies at the same rate as the disk particles', and large 
eccentricities and inclinations can get excited at these resonances. As planets sculpt the disk and cause its 
surface density to evolve over time, the location and strength of secular resonances will shift. Although 
radial drifts in the location of secular resonances might alter the details of how a planet depletes the disk 
as it excites particles into crossing orbits, this issue is likely of lesser importance when compared to the 
mean-motion resonances. 

In particular, a planet that is embedded in a self-gravitating disk can launch spiral density waves at 
its mean-motion resonances. Numerous examples of this phenomenon exist in Saturn's rings which exhibit 
density waves driven by orbiting satellites. The gaps in these rings reveal a history of angular momentum 
exchange between ring material and satellites, and similar exchanges are expected of planet-forming 
systems. When Neptune launches density waves at an exterior j + 1 '■ j mean-motion resonance, the disk 



exerts the torque Tj on the planet which opposes its radial migration ( Goldreich fc Tremaine 1980 ) 



^' = VdDjd-r " ~^^(-^fp^^^^e(a,n,) . (6) 

Here, ip — —8jfJ'p{rn)'^/5 is the planet's forcing function, fXp is the planet's mass in solar units, O is 
the disk's mean motion, D is the frequency difference from exact resonance which has the gradient 
rdD/dr ~ — 3(j + l)rj^, and the preceding quantities are evaluated at resonance r = (1 + j'^^Y^^cip 



where ap and Q.p are the planet's semimajor axis and mean motion (Hahn, Ward, & Rettig 1995). The 
dimensionless 'disk mass' is /id = ■na{ap)ai/MQ ~ 1.1 x \Q~^{Mo/Mi^) for the disk simulations that have 
cr(ap) — 1.3 X 10^"* Md/AU^. Figure H sums the torques Tj due to Neptune's outer j = 1 to 5 resonances 
for comparison with the non-resonant torque Tq that drives Neptune's orbit expansion. 

Evidently, the resonant disk torques can inhibit Neptune's outward migration if the disk admits 
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a wave response at the j ^ 3 resonances. It should be noted that these torques are also operative in 
non-self-gravitating disks (Lissauer fc Espresate 1998|), and similar torques also slow orbit expansion when 



resonance trapping is effective (Section 3.3). But when trapping is not of concern, the resonant torques in a 



non-gravitating disk shut off once particle eccentricities get excited ( Lissauer fc Espresate 1998 ). However 
in a self-gravitating particle disk, density waves transport the planet's forced disturbance downstream 
of the resonance in the direction of the planet's orbit. As long as wave action is sustained, particles at 
resonance maintain low eccentricities and the resonant disk torque can oppose planet-migration. 

There are two ways in which the propagation of density waves might be inhibited in particle disks. In 
a disk that is populated by comet-sized (or larger) planetesimals, dissipative forces such as gas drag or 



viscosity due to inter-particle collisions are insufficient to damp out density waves (e.g., Hahn, Ward, fc 



Rettig 1995). In this case, density waves reflect at a Q~barrier in the disk and return to the launch zone 



where they are reabsorbed by the particles at resonance (Toomrc 1969). This absorption of the returning 
waves' energy will steadily heat the disk and can eventually shut off subsequent wave generation. A second 
way to defeat waves is via stirring by larger bodies which also heats the disk and inhibits wave propagation. 
Consequently, the resonant torque that the particle disk exerts on Neptune will delay the onset of orbit 
expansion until the disk becomes too stirred to sustain density waves at its j ^ 3 resonances. 



7. Conclusions 

The existence of the Oort Comet Cloud as well as the Kuiper Belt suggest that there was once 
a residual planetesimal disk of mass ~ 10-100 M® in the vicinity of the giant planets following their 
formation. Further, any inefficiencies in the formation of the giant planets' cores implies additional disk 
mass. The eventual clearing of this planetesimal population involves a substantial exchange of orbital 
energy and angular momentum with the planets, implying that the present locations of the giant planets 
are not necessarily their formative ones. We have numerically simulated the evolution of a system of four 
giant planets embedded in a planetesimal disk of mass ranging from lOM® to 2OOM0. Our numerical 
simulations show a gradual increase in the mutual separation of the planets' orbits as the disk is dispersed 
via gravitational scattering by the planets. Higher disk masses yield planetary orbital migration that is 
faster and larger in magnitude. If planet-migration and resonance-trapping is invoked to explain the 
eccentricities of Pluto and its cohort of Kuiper Belt Objects at Neptune's 3:2 mean-motion resonance, then 
these simulations show that a disk mass of order Mo ~ 50 M^ is required to expand Neptune's orbit the 
requisite distance of Aa ~ 7 AU to pump up Plutino eccentricities to e ~ 0.3. Such an episode of planet 
migration implies that the initial Solar System was more compact in the past, with the Jupiter-Neptune 
separation having been smaller by about 30%. This finding also confirms the disk mass estimate previously 



obtained by Malhotra (1997) 



Our model disk-planet systems behave similarly to other disk systems that experience a gravitational 



or viscous torque (e.g., Lynden-Bell fc Kalnajs 1972, Lynden-Bell & Pringle 1974), which causes angular 



momentum to be carried radially outwards (in this application, by the outer three planets) while disk 
particles deliver mass radially inwards. However these particles tend to get ejected upon reaching Jupiter's 
orbit, which accounts for that planet's slight orbital decay. 

Since our simulations neglected the disk's self-gravity, collective effects such as density waves are 
precluded. A planet embedded in a self-gravitating planetesimal disk will tend to launch spiral density 
waves at its resonances. The torque due to wave generation is sufficient to oppose Neptune's orbit expansion 
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as long as the disk remains dynamically cold enough to admit a wave response from its j ~ 3 or higher 
resonances. Such an episode of wave generation will delay the onset of planet migration until the disk's 
wave response is defeated. 

The bulk of the disk particles deposited in the Oort Cloud originate in the vicinity of the Saturn- 
Neptune region of the solar nebula. Assuming that galactic tides and passing stars decouple particles from 
the planetary system when they achieve a semimajor axis of a > 3000 AU, and that these perturbations also 



remove about two-thirds of the Oort Cloud over the age of the Solar System ( Duncan, Quinn, fc Trcmaine 



1987 ), we estimate that about Mqc ^12 Mq of the Me, = 50 M^ disk is initially emplaced in the Oort 



Cloud, of which ^ 4 M^ will persist to the present age of the Solar System. 

Due to the fact that the disks simulated here were sparsely populated by particles having masses 
m = 0.01-0.2 Mg, their vigorous scattering caused the planets' orbits to evolve non-adiabatically such 
that resonance trapping of KBOs was inhibited. However previous studies have shown that adiabatic 
orbit expansion by Neptune can account for the abundance of eccentric KBOs that are known to orbit 
at Neptune's 4:3, 3:2, 5:3, and 2:1 mean-motion resonances (Malhotra 1995|). Unless there are (or were) 



additional unseen, distant perturbers, any primordial KBOs beyond Neptune's 2:1 resonance should reside 
in nearly circular, low-inclination orbits. 

The planet-migration/resonance trapping phenomenon might also have applications in extrasolar 
planetary systems. The most visible component of an extrasolar planetary system is likely its dust, which 
should be most abundant when planets and planetesimals are colliding, accreting, and eroding. Dusty 
circumstellar disks and rings are known to orbit the stars (3 Pictoris, Formalhaut, HR 4796A, 55 Cancri, 
and e Eridani (^mith fc Terrile 19"8^, [Greaves et al. 199^ , tTrilling fc Brown 1998| , [Koerner et al. 1998| , 



Holland et al. 1998 , Smith et al. 1998 ). In some of these systems, collisions and/or radiation forces will 
remove the observed dust on a timescale shorter than the age of the parent star. The presence of dust 
thus suggests an additional source — perhaps dust generation due to collisions by unseen planetesimals that 
reside in the disk. Any planets that might form within this environment will deplete the disk region that 
lies within their gravitational reach, which could account for these disks' central gaps. However an episode 
of planetesimal disk-clearing would also drive planet-migration, which can concentrate planetesimals at 
the outermost planet's exterior mean-motion resonances. Since the collision frequency and hence the dust 
generation rate varies as the square of the planetesimal density, one might speculate that this mechanism is 
also responsible for the formation of dust rings observed around e Eridani, Formalhaut, and HR 4796A. 
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Table 1: 

TABLE 1 

Oort Cloud Masses 
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Fig. 1. — Fractional variations in a particle's Jacobi integral J as a function of a particle's perihelion distance 
q averaged over 10^ years. Jupiter lies on a circular orbit at 5 AU and the integrator step size is At — 0.4 
years. Squares indicate particles having initial eccentricities of e — 0.1, inclinations i = 3°, and semimajor 
axes 4 < a < 50 AU, while crosses and circles are for eccentric particles having initial perihelia 0.05 < g < 6 
AU, a — & AU, and < z < 10°. Filled circles and squares indicate particles that approached closer than a 
Hill radius, or 0.35 AU, of Jupiter. 

Fig. 2. — A system consisting of Neptune on a circular orbit at a = 30 AU and 100 massless particles having 
initial semimajor axes 28.5 < a < 31.5 AU, eccentricities e = 0.05, inclinations i = 0.025 radians, and 
Jacobi integrals 2.997 < J < 2.999 is evolved for 5 x 10^ years. The solid curve gives the swarm's average J 
with vertical bars indicating the swarm's standard deviation aj. The characteristic timescale to diffuse into 
Uranus-crossing orbits having J < Ju that are dynamically forbidden is td ~ 2.5 x 10^ years. 

Fig. 3. — (a) The fractional error in the relative velocities 8V/V versus periapse q (in units of Neptune's 
Hill radius Rh) of 500 particles after scattering off Neptune. The particles' initial heliocentric orbits were 
a ~ 30 AU, e = 0.05, and i = 0.025 radians. Particles that trigger the two-body close encounter algorithm 
are indicated by a dot. (b) A histogram of the z-component of the specific angular momentum errors 5(.p for 
each scattering event. The total specific angular momentum exchanged between the planet and the particle 
swarm is L = | ^-^pj = 1.27 AU^/yr which has an rms sum (X]^p)^^^ = 0-64 AU^/yr and an rms error 
5L = 0.020 AUVyr . 

Fig. 4. — The semimajor axes of the the giant planets while embedded in planetesimal disks of mass Md = 10, 
50, 100, and 200 M0. The boundaries of the grey regions denote the planets' perihelia and aphelia distances. 

Fig. 5. — The eccentricities e versus semimajor axes a at logarithmic time intervals for the M^i = 50 Mq 
system. Small dots indicate scattered particles that have passed within a Hill radius of a planet and crosses 
indicate particles that have not encountered a planet. Large dots denote the planets and the vertical dashes 
indicate the location of Neptune's outer four mean-motion resonances. Orbits lying above the left curve 
have perihelia inside of Saturn's orbit and those above the right curve have perihelia interior to Neptune. 

Fig. 6. — The simulation of the AIo = 50 M0 system extended out to i = 5 x 10^ years. Grey indicates 
perihelia and aphelia distances. Note that Uranus and Neptune pass through a 2:1 mean-motion resonance 
at i = 3.05 X 10^ years, which results in brief eccentricity excitation. 

Fig. 7. — The radial displacement Aa versus disk mass Md for each planet after t — ?> x 10^ years. 

Fig. 8. — The solid curve is the time-averaged torque Tq on Neptune, in units of that planet's current angular 
momentum/orbital period ratio L/P, and is plotted versus the disk mass Me,- The dotted curves indicate 
sums of the resonant disk torque contributions, Eq. y. 

Fig. 9. — The fraction of the active disk that is ejected into hyperbolic orbits, //i, versus time t for the 
Mu = 100 M0 simulation. The disk fraction that is scattered into the Oort Cloud with a > 3000 AU is fsk, 
and /g<3.5 is the instantaneous disk fraction having perihelia q < 3.5 AU. Note the near constancy of the 
ratio fzk/fh — 0.4. Curves for the other runs are quite similar. 
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Fig. 10. — The extrapolated mass of the Oort Cloud mass Mqc as a function of the initial disk mass Mr,. 
The dashed curve includes contributions by massless test particles originating in the 4 < a < 10 AU part of 
the disk that is modeled separately. However each one of these particles, when deposited in the Oort Cloud, 
are assumed to contribute the same individual masses as their compatriot particles that are employed in the 
simulations of Fig. ^. 



Fig. 11. — The local disk fraction deposited in the Oort Cloud versus initial semimajor axis a for the Md = 50 
Mq simulation. The data are obtained at time i = 3 x 10^ years and then extrapolated to a fully evolved 
state by multiplying by f'^k/ fzk — 1-97. The grey bars show the extent of radial migration by each planet. 
The 4 < a < 10 AU component (dashed curves) is obtained from a separate integration, and N^''^ errors are 
assumed. 



Fig. 12. — Eccentricity e versus semimajor axis a at time i = 5 x 10^ years for the Mo = 50 M^ simulation. 
Large dots indicate Uranus and Neptune, small dots indicate scattered particles that have passed within a 
Hill distance of a planet, and crosses denote particles that have not had a close planetary encounter. The 
dashed lines indicate Neptune's four outermost mean-motion resonances, and the ej=3 curve satisfies Eq. 
(0) with i = 0. Boxes denote observed KBOs having well-determined orbits (from Marsden (1998)), and 
scattered object 1996 TLgg is indicated. 
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